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ABSTRACT 


STADAN Minitrack network data can be adequately an- 
alyzed using an optically determined '’standard orbit” of 
5-1/4 days based primarily upon GEOS I flashing- lamp se- 
quence observations recorded by STADAN and SPEOPT 
cameras. A technique has been devised to uncouple the 
refraction-dependent from the refraction-independent ef- 
fects in the Minitrack residuals. Results indicate the 
Minitrack system as a whole to be free of any strong bias. 
Station-dependent systematic errors were not removed 
by the usual calibration procedures, and accordingly an 
optically- determined ’’standard orbit” would be a preferable 
method of calibrating individual Minitrack antenna arrays. 
A Minitrack orbit for the same time period as the optical 
orbit was obtained, and position differences between the two 
orbits produced an rms deviation of approximately 165 
meters— a favorably small figure. 



SUMMARY 


STADAN Minitrack network data are analyzed using an optically determined 
"standard orbit" 5-1/4 days in length as a reference orbit. This "standard orbit" 
has an rms fit of 3"0. It is shown that this is adequate for evaluating Minitrack 
data. The "standard orbit" is based primarily upon GEOS I flashing-lamp sequence 
observations recorded by STADAN and SPEOPT MOTS 24" and 40", SAO Baker- 
Nunn, and USAF PC- 1000 cameras. 

A technique is devised to uncouple the refraction- dependent from the 
refraction-independent effects in the Minitrack residuals. The mean value of the 
refraction-dependent residuals exhibited a shift of -0.2 x 10 ' 3 in the region of 10° 
to 30° elevation. The shift is insignificant at higher elevations. The mean value of 
the refraction-independent residuals is approximately 0.02 x 10 ~ 3 with a standard 
deviation of 0.2 x 10" 3 . This small mean value indicates that the Minitrack system 
when viewed as a whole appears relatively free of any strong bias. The Minitrack 
data were not included in the reference orbital solution. Analysis of the residuals 
by individual stations indicates that there may be some station -dependent system- 
atic errors present. Since these systematic errors were not removed by the usual 
calibration procedures, it is suggested that the use of an optically- determined 
"standard orbit" would be a preferable method of calibrating individual Minitrack 
antenna arrays. 

A Minitrack orbit for the same time period as the optical orbit was obtained 
with an rms fit of 0.19 x 10 ~ 3 . The Minitrack data used in this particular solution 
were not corrected for ionospheric or tropospheric refraction and were not 
weighted for elevation effects. Data down to 15° elevation were used in the solution. 
Position differences between the Minitrack and optical orbits produced an rms de- 
viation of approximately 165 meters. This low figure is due in part to the use of an 
accurate gravity model and of station positions consistent with this gravity model. 


ABBREVIATIONS USED IN THIS DOCUMENT 


MOTS Minitrack Optical Tracking System 

SAO Smithsonian Astrophysical Observatory 

SPEOPT Special Optical Tracking System 

(Same cameras used by MOTS) 

STADAN Satellite Tracking and Data Acquisition Network 


TRACKING- STATION CODE NAMES 

A detailed description of the names, general locations, and specific coordinates of the tracking 
stations corresponding to the adopted six-character station code names used in the various tables 
and figures of this document is given in Appendix D. 
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INTERCOMPARISON OF THE MINITRACK 
AND OPTICAL TRACKING NETWORKS USING 
GEOS I LONG-ARC ORBITAL SOLUTIONS 


by 

James G. Marsh and Clarence E. Doll, Jr. 
Goddard Space Flight Center 

and 

Richard J. Sandifer and Wayne A. Taylor 
Wolf Research and Development Corporation 


INTRODUCTION 

The primary objectives of the GEOS I project are to use observations of GEOS I for the pur- 
poses of connecting geodetic datums and defining the earth 1 s gravitational field to a specified de- 
gree of accuracy, improving positional accuracies of satellite tracking sites, and evaluation and 
calibration of tracking equipment (Reference 1, pp. 1, 6). This report is concerned with the last- 
named aspect of evaluation and calibration as applied to the Minitrack system. 

One particular phase of the Minitrack/optical data intercomparison is presented in this report. 
GEOS I long-arc orbital solutions determined by optical observations only are used as a "standard 
orbit/’ The Minitrack residuals used in this analysis are calculated on the basis of this optically- 
determined T ’standard orbit.” As a result, the Minitrack residuals are calculated on the basis of an 
adjusted orbit which is not influenced by the Minitrack data set itself. Assuming that there are no 
significant systematic differences between the two data types, or— perhaps equally important— no 
significant systematic differences between the tracking networks which recorded the two data types, 
this type of inter comparison can provide a useful tool for analyzing and evaluating the Minitrack 
data and network. It will be shown in a later section that the geometry of the reception pattern of 
the Minitrack antenna arrays provides a possible method of uncoupling refraction- dependent effects 
from refraction-independent effects contained in the Minitrack residuals. In order to analyze the 
data from this viewpoint, the Minitrack residuals have been calculated from data which have not been 
preprocessed. As a result, the full effects of both tropospheric and ionospheric refraction are in- 
cluded in the residuals. By analyzing the data in this manner, it may be possible to determine how 
successfully the uncoupling process can be performed. If the process is successful, there is the 
added advantage that the data do not contain erroneous refraction effects which may arise from over- 
correcting or under-correcting the data for refraction. This problem is particularly sensitive be- 
cause of the large uncertainties contained in present knowledge of the state of the ionosphere. The 
analysis in this report wilhbe restricted to Minitrack data with all refraction effects included. 
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The last section herein before T 'conclusions’ T presents an inter comparison of two 5-1/4-day 
trajectories. One trajectory was computed on the basis of initial conditions adjusted by optical 
data only. The second trajectory was computed on the basis of initial conditions adjusted by Mini- 
track data only. 

Using the optic ally -determined orbit as the standard, the errors in position as determined by 
the Minitrack data had a standard deviation of approximately 165 meters. The Minitrack solution 
utilized data down to approximately 15° elevation. These data were not corrected for tropospheric 
or ionospheric refraction, and were not weighted to compensate for any deterioration of the data at 
the lower elevations. The reduction of the standard deviation to the relatively low value of 165 meters 
may be attributed in part to the use of a sophisticated gravity model (SAOM-1 model— Reference 2) 
and the transformation of all tracking site locations to a datum (SAO C- 5— Reference 2) using the 
center of mass of the earth as its origin and consistent with the gravity model used. 

FEASIBILITY OF USING OPTICALLY-DETERMINED ''STANDARD ORBITS” 

TO EVALUATE AND CALIBRATE MINITRACK DATA 

The feasibility of using optical observations of artificial satellites for calibration purposes 
has been previously noted by V. R. Simas.* The large volume of high-precision optical observa- 
tions of GEOS I currently available offers the first opportunity to analyze in depth the data obtained 
from the Minitrack network of the same satellite and to explore the feasibility of utilizing such 
optical data for calibration purposes. As noted in Appendix A, the upper precision limit obtainable 
at the most basic level of the electronics of the Minitrack system (the phase-measurement level) 
is about 0.36 electrical degrees. This is equivalent to a precision of approximately 5” in terms of 
ability to measure angular arc in space under optimum conditions at the zenith. A more realistic 
figure for the overall expected precision of the Minitrack system at the equipment level is probably 
around 20” in terms of angular direction in space (Reference 3, p. 84). The precision of the SAO 
optical observations is estimated to be approximately 2” (Reference 2, p. 43) for a single observa- 
tion. Observations from the MOTS, PC- 1000, and other cameras used in the GEOS I program have 
an inherent precision which is slightly better than the SAO figure. Whether or not this inherently 
greater accuracy is actually realized will depend to a great extent upon the care exercised in re- 
ducing and preprocessing the data. Taking the SAO figure of 2” as an approximate gauge of the 
overall accuracy of an individual optical observation, there is sufficient justification for using op- 
tical observations as a standard to be used in analyzing and evaluating Minitrack network data pro- 
vided that there exist no significant systematic differences between the two data types. In a later 
section, this will be shown to be the case when considering the Minitrack network as a whole, al- 
though there appear to be systematic problems at individual Minitrack sites. Even if there exist 
significant systematic differences between the two data types, there still is a strong justification 
for using optical data as a ’’standard,” since all detectable error sources in this data type as exist 
in star catalogues, plate measuring equipment, optical refraction, etc., have had the advantage of 
decades of thorough study by competent investigators, so that the error limits of these sources are 
well known. 


Unpublished notes. Sections 3.4.4 and 3-4.5, p- 54. 


A certain amount of caution should be exercised when considering the MOTS optical observations 
as used in the '’standard orbit,” since the MOTS cameras are used also for calibrating the Minitrack 
system and any undetected systematic errors in the MOTS instrumentation will probably be propa- 
gated into the Minitrack calibration coefficients. However, since the MOTS cameras are co-located 
with the center of the Minitrack antenna ’’fine beam” array, these optical observations may prove 
useful in verifying suspected systematic errors which have been detected by the Minitrack residuals. 

One particularly useful conclusion of the current investigation is that an optically- determined 
’’standard orbit” has great value in detecting systematic errors at individual Minitrack sites. K 
these individual systematic errors can be identified and removed, the overall result will be a cor- 
responding reduction in the rms of fit of orbits as determined by Minitrack data only. It is sug- 
gested also that the use of an optically -determined ’’standard orbit” would be a preferable method 
of calibrating the individual Minitrack sites. The current method of calibration by means of simul- 
taneously photographing an airborne flashing light with the MOTS cameras and detecting the 136.5- 
Mc radio signal with the Minitrack antenna arrays is not capable of detecting all systematic errors. 
Some obvious examples are errors in site location and timing errors. At best, the present cali- 
bration method can account only for local systematic errors. 

SUMMARY OF OPTICAL DATA USED IN THE "STANDARD ORBIT" SOLUTION 

The GEOS I optical and Minitrack data used in the analysis described in this report were ob- 
tained from the NASA Space Science Data Center located at the Goddard Space Flight Center. 

The ’’standard orbit” used in the calculation of the Minitrack residuals in the current analysis 
was determined completely from optical observations of GEOS I. Minitrack observations were not 
used to adjust the standard orbit which was used as a basis for the calculation and subsequent anal- 
ysis of the Minitrack residuals. 

For the purpose of calibrating electronic equipment, the ’’active” GEOS I optical data are 
particularly useful as a high-precision data set. In addition to the reasons noted earlier, other 
factors contributing to this precision are: 

1. The use of a stable on-board clock to trigger the optical beacon flash sequences permitted 
the determination of the time of observation to millisecond accuracy. 

2. The short duration (approximately 1.3 milliseconds) of the optical beacon flash sequences 
permitted the tracking cameras to record the observations as point images rather than as 
a streak against a background of reference stars, thus making it possible for the right as- 
censions and declinations of the flashes to be determined to higher precision than normally 
obtained. 

In general, the reduction methods associated with using reference stars as a means of deter- 
mining angular position are the most accurate available among all current tracking systems. 

The epoch of the standard orbit was 01 H 38 M 225000 1^3 1/1965 UTC time. The epoch was chosen 
at the beginning of the selected optical data set which extended to 06 H 45 M 1/5/1966 UTC, a total arc 
length of approximately 5-1/4 days covering 63 orbital revolutions of the spacecraft. The resulting rms 
of fit of this particular standard orbit was 3 ’.’12 based upon 1059 observations summarized in Table 1. 
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Table 1 


Summary of GEOS I Optical Observations Used in the Orbital Solution 
for the Period 1 H 38 M 12/31/65 or 6 H 45 M 1/5/66 (63 Orbital Revolutions). 





No. of Observations 

Network 

Station 

Camera Type 

a 


8 

a and 8 

Type 

No. Passes/ 
No. Flash Seq. 

SAO 

1 ORGAN 

Baker-Nunn 

1 


1 

2 

Passive 

i/ 


1MAUIO 

Baker-Nunn 

1 


1 

2 

Passive 

i/ 


1NATOL 

Baker-Nunn 

4 


4 

8 

Passive 

4/ 


OSLONR 

Baker-Nunn 

2 


2 

4 

Passive 

2/ 


AUSBAK 

Baker-Nunn 

2 


2 

4 

Passive 

2/ 


1SHRAZ 

Baker-Nunn 

1 


1 

2 

Passive 

2/ 


1SPAIN 

Baker-Nunn 

3 


3 

6 

Passive 

3/ 


1TOKYO 

Baker-Nunn 

6 


6 

12 

Passive 

6/ 


1VILDO 

Baker-Nunn 

1 


1 

2 

Passive 

1/ 


1JUPTR 

Baker-Nunn 

14 


14 

28 

Active 

1/2 


AGASSI 

K-50 

5 


5 

10 

Active 

i/i 


TOTAL 


40 + 40 = 80 

(21 +21 =42 total passive) 



SPEOPT 

1COLBA 

MOTS 40" 

81 


83 

164 

Active 

8/13 


1JUM40 

MOTS 40" 

7 


15 

22 

Active 

3/3 


1BERMD 

MOTS 40" 

39 


45 

84 

Active 

5/7 


1PURIO 

MOTS 40" 

7 


7 

14 

Active 

i/i 


1DENVR 

MOTS 40" 

28 


42 

70 

Active 

4/6 


1JUM24 

MOTS 40” 

10 


16 

26 

Active 

3/3 


TOTAL 


172 

+ 

208 

= 380 

Active 


STADAN 

1FTMYR 

MOTS 40" 

41 


41 

82 

Active 

4/6 


1BPOIN 

MOTS 40” 

23 


30 

53 

Active 

5/6 


1GFORK 

MOTS 40" 

11 


15 

26 

Active 

3/3 


1MOJAV 

MOTS 40" 

12 


13 

25 

Active 

2/2 

[ 


TOTAL 


87 

+ 

99 

= 186 



USAF 

HUNTER 

PC-1000 

30 


29 

59 

Active 

5/5 


SWANIS 

PC-1000 

7 


7 

14 

Active 

1/1 


GRDTRK 

PC-1000 

0 


7 

7 

Active 

1/1 


ANTIGA 

PC-1000 | 

13 


13 

26 

Active 

2/2 


SEMMES 

PC-1000 

30 


30 

60 

Active 

4/5 


CURACO 

PC-1000 

20 


20 

40 

Active 

3/3 


HOMEST 

PC-1000 

47 


47 

94 

Active 

4/6 


JUPRAF 

PC-1000 

12 


5 

17 

Active 

2/2 


BEDFRD 

PC-1000 

11 


11 

22 

Active 

2/2 


ABERDN 

PC-1000 

37 


37 

74 

Active 

6/6 


TOTAL 


207 ■ 

+ 

206 

= 413 


21/ Passive 
70/86 Active 


Total of All Observations = 1059 
Total Passive Observations = 42 


91 Total Station-Passes 
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A more detailed description of the names, general locations, and specific coordinates of the track- 
ing stations corresponding to the adopted six- character station code names listed in Table 1 and 
other tables and figures of this report is given in Appendix D. The position and velocity vectors 
for the epoch 01 H 38 M 22^000 12/31/1965 UTC which were adjusted on the basis of the 1059 optical 
observations are given in Table 2. A summary of optical data coverage by time is shown in 
Figure 1. 


Table 2 


Final Adjusted Position and Velocity Vectors at 
Epoch for the Optical ’’Standard Reference Orbit.” 


Optically Adjusted 
Position Vector 

X: +5,690,537.7 meters 
Y: +1,474,538.5 meters 
Z: +6,013,442.9 meters 

Epoch: 

No. of Observations: 
Arc Length: 

Rms of Fit: 


Optically Adjusted 
Velocity Vector 


X: -4,685.6198 meters/sec. 
Y: +3,849.4695 meters/sec. 
Z: +2,939.1210 meters/sec. 

01 H 38 m 225000 12/31/1965 UTC 

1059 

5-1/4 days 
3712 


For comparable arc lengths, this rms of 

OPTICAL 

DEC 

- + - 

31 

■ 

1 ■ 

. J 

JAN 

1 

■ 

1 

fit was not of the highest quality. A more typ- 

MINITRACK 

ical figure for GEOS I for this arc length would 

mm T m 

l-Ll 1111 LLl-L 

-LLULL. 

■ ■1 

Mini 


M ■ 

II 1 1 1 

JJ..I 1 L 

■ ■■ 

JLL-1 1 1 


be approximately 21’0. However, this particu- 
lar arc had been used for intercomparison with 
other instrumentation types and it was felt that 
a multiple-instrument intercomparison would 
prove useful. In addition, the 3 T . T 12 rms of fit 
still places the optical data set well within 
the precision limits required for a Minitrack 
intercomparison. 


JAN 2 


JAN 3 


OPTICAL 

MINITRACK 


OPTICAL 

MINITRACK 


■ 

■ 

1 

■J 

■ ■ 


■ ■ 

m^mmm ■ ■■ 

tgzm 

■ ■ ■ 

1 


mu juii mu mjii 


1 1 l-U. 

Mill- 

j-u 1 1 


JAN 4 


JAN 5 


t- 

■ 1 

* * 

■ 

i 

■ 


L. 

II 


ll> 

■ 


Him 

mu 


liJLLL 

1 1 1 U 

mn. 


12 


18 


12 


18 


24 


Figure 1— Summary of GEOS I optical and Minitrack 
coverage for the period 1 H 38 M 12/31/65 to6 H 45 M 1/5/66 
UTC (63 orbital revolutions). 


The 1059 optical observations are sum- 
marized in Table 1 by network, station, and 
type, i.e., whether the data were active or pas- 
sive. In addition, a further breakdown is given 
of the active observations according to the 

number of passes over a particular station and the total number of flash sequences (maximum of 
7 right ascensions and 7 declinations per flash sequence) per station. These data are more 
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descriptive of the available station coverage and geometrical strength if they are described in terms 
of "station-passes." 

The final tally reduces the total of 1059 optical observations to 86 active flash sequences dur- 
ing 70 "active station- passes" over the various optical stations and 21 "passive station-passes." 

The total number of "station-passes" was 91. The 21 "passive station-passes" were all obtained 
from the SAO network. Of the 70 "active station-passes," only 2 were obtained from the SAO net- 
work. The remaining 68 "active station-passes" were distributed as follows: 

30 passes USAF Network PC- 1000 cameras 
24 passes SPEOPT Network MOTS 40" cameras 
14 passes STADAN Network MOTS 40" cameras 

The 68 passes from the Air Force, SPEOPT, and STADAN networks were all located on the 
North American continent or its near vicinity (Figure 2). In addition, the only two "active station- 
passes" from SAO, at Agassiz, Massachusetts, and Jupiter, Florida, are on the North American 
continent (Figure 3). At least one station-pass was obtained from every SAO station shown in Fig- 
ure 3 except Olifantsfontein in South Africa. The rms of fit of 3712 for the 5-1/4 days of optical 
observations was obtained using the SAO M-l gravity model modified by the GEOS I resonant 


132 ° 120 ° 108 ° 96 ° 84 ° 72 ° 60 ° 54 ° 



Figure 2 — Location of USAF, SPEOPT, and STADAN camera sites 
whose observations were used in the orbital solution. 


6 




harmonics (c l3 12 , C 14 12 , C ls 12 , s l3 12 , s 14 12 , S is 12 ) as redetermined by SAO (Reference 2, p. 2; 
Reference 4, p. 5; and Appendix C). All tracking station locations were transformed to the SAO 


C-5 system (Reference 2, p. 2; Appendix D). 


SUMMARY OF MINITRACK DATA AVAILABLE FROM THE SPACE SCIENCE DATA 
CENTER FOR THE PERIOD OF INTEREST, AND DATA REJECTED BEFORE ANALYSIS 

The Minitrack data as obtained from the NASA Space Science Data Center consist of 1 to 3 
pairs of direction cosines per station-pass. These data in turn have been reduced by Goddard from 
approximately 30 pairs of phase -difference measurements per pass as transmitted to Goddard from 
the various Minitrack stations. Goddard transforms these data to direction cosines and reduces 
them to approximately 1 to 3 data pairs of direction cosines by means of a polynomial smoothing 
process. This is the only data preprocessing performed before the data are transmitted to the 
Space Science Data Center or used by Goddard T s own orbit differential correction program. The 
rms of fit of the polynomial smoothing process is slightly less than 0.1 x 10" 3 on the average.* 

During normal operations, the Minitrack stations do not track below 20° zenith distance except 
during early launch phases. In the case of GEOS I, an exception was made and tracking was performed 
almost to the horizon. For the current analysis, Mini track observations are available down to 15° 
elevation, so that there is a sufficient span of data in that dimension to investigate refraction effects. 

The Minitrack data originally available for analysis from the Space Science Data Center for the 
period of interest consisted of 482 observations, summarized by station in the first column of Table 3. 


*Watkins, E. R., Jr., private communication. Operational Computing Branch, Computation Division, Goddard Space Flight Center. 
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Table 3 


Summary of Minitrack Data Rejected for Analysis Purposes for the 
Period 1 H 38 M 12/31/65 to 6 H 45 M 1/5/66 (Approximately 63 Orbital Revolutions). 




{ 




m 


Station 

Total No. 

of Obs. 
Available 

No. of Obs. 
with Residuals 
Greater 

Than 0.5 x 10" 3 

Total No. 
of Obs. 
Rejected 

% of Obs. 
with Residuals 
Greater 

Than 0.5 x 10“ 3 

Total No. 

of Obs. 
Available 

No. of Obs. 
with Residuals 
Greater 

Than 0.5 x 10“ 3 

Total No. 
of Obs. 
Rejected 

% of Obs. 
with Residuals 
Greater 

Than 0.5 x 10“ 3 

1. BPOINT 

50 

1 

4 

2.0 

50 

2 

4 

4.0 

2. COLEGE 

25 

3 

4 

12.0 

25 

2 

4 

8.0 

3. FTMYRS 

21 

2 

4 

9.5 

21 

3 

4 

14.3 

4. GFORKS 

29 

2 

2 

6.9 

29 

1 

2 

3.4 

5« JOBURG 

2 

0 

0 

0.0 

2 

0 

0 

0.0 

6. LIMAPU 

6 

0 

0 

0.0 

6 

0 

0 

0.0 

7. MOJAVE 

38 

1 

2 

2.6 

38 

1 

2 

2.6 

8. NEWFLD 

38 

0 

2 

0.0 

38 

1 

2 

2.6 

9. OOMERA 

6 

1 

2 

16.7 

6 

1 

2 

16.7 

10. QUITOE 

4 

0 

0 

0.0 

4 

0 

0 

0.0 

11. SNTAGO 

4 

0 

0 

0.0 

4 

0 

0 

0.0 

12. WNKFLD 

18 

3 

3 

16.7 

18 

2 

3 

11.1 

i 

241 

13 

23 

5.4% 

241 

13 

23 

5.4% 


m 

241 

13 

23 

5.4% 

TOTAL 

482 

26* 

46 

5.4% 


5.4% of data. 



After a rejection criterion (described in the next paragraph) was applied to the data, 436 data points 
remained for the analysis. These data are summarized in Table 4. The 436 data points comprise 
158 station-passes, compared to the 91 station-passes of optical data (Table 1). For the Minitrack 
data also, the geometry is strongly oriented toward the North American continent, with 29 passes 
at Blossom Point, Maryland, 26 passes at St. John's, Newfoundland, and 25 passes at Mojave, 
California. At least one observation pair was obtained from every Minitrack station (Figure 4) 
except the one located at Tananarive. 

A cutoff criterion of 0.5 x 10 ~ 3 in direction cosine residuals was used in the rejection of data 
to be used in the analysis. This particular figure was chosen for convenience in the analysis of the 
residuals and to correspond to approximately 2-1/2 times the rms of the orbital fit. For data 


Table 4 

Summary of Minitraek Observations Used in Intercomparisons for the 
Period 1 H 38 M 12/31/65 to 6 H 45 M 1/5/66 (Approximately 63 Orbital Revolutions). 




No. of Observations 



Station 

l 

m 

No. of 

Location 




Passes 


1 . 

BPOINT 

46 

46 

29 

Blossom Point, Maryland 

2. 

COLEGE 

21 

21 

18 

College, Alaska 

3. 

FTMYRS 

17 

17 

12 

Fort Myers, Florida 

4. 

GFORKS 

27 

27 

19 

East Grand Forks, Minnesota 

5. 

JOBURG 

2 

2 

2 

Johannesburg, South Africa 

6. 

LIMAPU 

6. 

6 

4 

Lima, Peru 

7. 

MOJAVE 

36 

36 

25 

Mojave, California 

8. 

NEWFLD 

36 

36 

26 

St. John*s, Newfoundland 

9. 

OOMERA 

4 

4 

4 

Woomera, Australia 

10. 

QUITOE 

4 

4 

3 

Quito, Ecuador 

11. 

SNTAGO 

4 

4 

4 

Santiago, Chile 

12. 

WNKFLD 

15 

15 

12 

Winkfield, England 





158 Total Station-Passes 


TOTAL 

218 + 218 

= 436 



Rms orbital fit using Minitrack data only - 0. 19 x 10 3. 

All Mini track data rejected whose residuals from optically-determined orbit exceeded 0.5 x 10 (5% of smoothed data available from 

data center). 
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Figure 4— STADAN 136-Mc Minitrack interferometer site locations. 


rejection purposes also, the "standard orbit" was used as the basis for calculating the Minitrack 
residuals. 

Approximately 5 percent of all Minitrack data available from the NASA Space Science Data 
Center during the period of analysis exceeded this figure. Thus the percentage of data rejected on 
the basis of the 0.5 x 10" 3 cutoff criterion is statistically consistent with the assumption of a normal 
distribution of the residuals (see residual plot in Figure 7). Minitrack data rejected on the basis 
of this criterion are summarized by station in Table 5. The magnitude of the residuals enclosed 
by rectangles in this table was less than the criterion, but they were omitted from the analysis be- 
cause their counterpart exceeded the criterion. One residual pair from Blossom Point, Maryland, 
at 17° elevation was erroneously rejected because of a faulty computer printout on the first pass 
through the data. A final summary of the total number of original observations and number of ob- 
servations rejected by station is shown in Table 3. The total number of observations rejected from 
the analysis was 46 of a total of 482, or 9-1/2 percent of the total data available. This is higher 
than the 5.4 percent of the data which actually exceeded the rejection criterion of 0.5 x 10~ 3 , because 
both data points of a pair were rejected even though only one of the two points may have exceeded 
the criterion. The rejection of 16.7 percent of the data from Woomera, Australia, is statistically insig- 
nificant because of the small number of observations available from that station (6 pairs). Of perhaps 
more significance is the lOpercent rejection rate for College, Alaska. It will be shown later that there 
is a strong bias in the residuals from this station which may account for the high rejection rate. 
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Table 5 


Table of Residuals of Mioitrack Data Rejected for Analysis Purposes 
(Rejection Criterion— 0.5 x 10 " 3 ). 


Station 
G FORKS 

BPOINT 

COLEGE 

FTMYRS 




2.54 

6.02 

507.12 


.19 


.03 


.29 


18.15 

17.97 

17.40 

.51 

.74 

.61 


.27 


.41 


x 10 3 
x 10 “ 3 

X 10“ 3 
X 10“ 3 
x 10“ 3 
X 10“ 3 

x 10" 3 
x 10" 3 
x 10" 3 
X 10“ 3 

x 10" 3 
x 10" 3 
x 10" 3 
x 10" 3 


MOJAVE 

17.22 

x 10 3 

1 -06 | 



1 - 12 1 

X 10“ 3 

.54 

NEWFLD 


.23 

X 10“ 3 

.52 



.12 

X 10“ 3 

.61 

WNKFLD 

132.50 

X 10“ 3 

113.17 


187.89 

x 10“ 3 

47.89 


1.99 

x 10“ 3 

rrn 

OOMERA 


.62 

x 10“ 3 

I .25 | 


1 

.38 

) X 10“ 3 

.62 


Am 


1.74 

DE] 

452.86 

.61 

rosi 

.83 


.30 


23.78 

17.56 


.15 

T09~ 

.55 

.58 

.59 


x 10“ 3 
x 10" 3 

x 10" 3 
x 10" 3 
x 10“ 3 
x 10" 3 

x 10“ 3 

x 10 -3 
x 10" 3 
x 10" 3 

x 10” 3 
x 10" 3 
X 10 " 3 
X 10" 3 

X 10 " 3 
x 10" 3 

X 10“ 3 
X 10 " 3 

x 10 " 3 

x 10-3 
X 10-3 


x 10 
x 10 


Elevation 


71° 

44 

-56 

23 

17 

21 

23 

20 

19 

28 

50 

23 

28 

27 

22 

34 

17 

70 

70 

61 

62 

47 

63 


Any residual enclosed by a rectangle (except the BPOINT residuals at 17 elevation) is less than the rejection criterion but was 
omitted from the analysis because its counterpart exceeded the criterion. Although both BPOINT residuals at 17 elevation were less 
than the rejection criterion, they were erroneously rejected. 


An orbit is described later (in the last section of this report before T ’Conclusions ") which is 
adjusted on the basis of Minitrack data only. These data were not corrected for tropospheric or 
ionospheric refraction effects and were not weighted for elevation effects in the orbital adjustment 
process. The resulting rms of fit was 0.19 x 10~ 3 . It is significant that when the Minitrack resi- 
duals were recalculated on the basis of the same orbital arc as adjusted by Minitrack data only, 
there were no additional residuals encountered falling outside the 0.5 x 10“ 3 criterion originally 
used to reject Minitrack residuals calculated on the basis of the "standard orbit." 


GEOMETRY OF THE MINITRACK "EQUATORIAL” 

AND "POLAR” MODES OF TRACKING 

The 136-Mc Minitrack interferometer tracking system and its tracking geometry are described 
in Appendix A. The tracking geometry is summarized in Figure 5. There are two possible modes 
of tracking, the "equatorial" tracking mode and the "polar" tracking mode. In the equatorial mode, 
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ZENITH 


The Mini track Observations are 


ZENITH 



( Long beam oriented N -S) 

At- Dependent upon long beam only 

Am - Dependent upon long beam & refraction 


t (- cos a) 
and m ( = cos/3) 


50 ° 



( Long beam oriented E -W) 

At - Dependent upon long beam & refraction 
Am - Dependent upon long beam only 


Figure 5— Approximate reception pattern of the fine beam of the 136-Mc 
Minitrack antenna array in the equatorial and polar tracking modes. 


the "fine beam" is oriented in a north-south direction, with the long portion of the beam stretching 
50° on each side of the zenith and 5° on each side of the zenith in the east- west direction relative 
to the station center C. When tracking in the polar tracking mode, the long portion of the "fine 
beam" is oriented east-west and the narrow portion oriented north-south. The particular mode of 
tracking is usually determined by the direction in which the satellite is approaching the station. It 
is possible to switch electronically from one mode to the other fairly quickly. Several passes of 
GEOS I during the period of analysis were tracked in both modes during the same pass. 

The observations available as data are the direction cosines l and m, which are equal respec- 
tively to the cosines of the angles a and fi as indicated in Figure 5. A vector directed toward the 
satellite is indicated in this figure in both the equatorial and polar tracking modes. Observational 
data are usually available only when the satellite is located inside this fan-shaped fine beam, al- 
though there are occasional data outside this pattern in the "side lobes." Observations made 
within the fine beam are considered to be more accurate than observations made in the side 
lobes. 

The method of obtaining calibration coefficients for the individual Minitrack stations is de- 
scribed in References 5 through 8. The method employed by Goddard Space Flight Center in trans- 
forming the raw Minitrack phase measurements to direction cosines, making use of these 
calibration coefficients and other relevant information, is described in References 9 through 12 
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and 13a, b, and c.* The corrections for ionospheric and tropospheric refraction as employed by 
the operational orbit determination program are described in References 14a and b and 15. 

Refraction corrections were not applied to the measurements used in the analysis performed 
in this report, as mentioned also in the preceding section. 

The Minitrack antenna arrays are laid out by a method of leveling corresponding to the astro- 
nomical zenith. The raw phase measurements are referenced to an "electrical" zenith defined by 
the antenna system. However, the calibration equations referred to above transform the raw phase 
measurements to direction cosines relative to the geodetic zenith. 

USE OF THE MINITRACK GEOMETRY TO UNCOUPLE REFRACTION-DEPENDENT 
FROM REFRACTION-INDEPENDENT EFFECTS IN THE MINITRACK RESIDUALS 

By considering the geometry of the antenna pattern in Figure 5 more closely, it is seen that, 
in the equatorial tracking mode, the angle a will never be less than 85° or more than 95° and the 
angle /3 will be approximately equal to the elevation angle or its supplement. In the polar tracking 
mode, the situation is reversed and fi will always be between 85° and 95° and a will be approxi- 
mately equal to the elevation angle or its supplement. For the purposes of the analysis in this re- 
port, it will be assumed on the basis of this configuration that any systematic errors due to re- 
fraction are not present in At and Am , , (At = residuals in the direction cosine l when 

tracking in the equatorial tracking mode; Am olar = residuals in the direction cosine m when track- 
ing in the polar tracking mode). Whatever systematic effects may be present in Af eq and Am polar 
will be called a "beam only" effect. It will be assumed that At , and Am , i.e., the residuals in 
the direction cosine l when tracking in the polar tracking mode and in the direction cosine m when 
tracking in the equatorial tracking mode, will contain both systematic "beam only" effects and sys- 
tematic "refraction" effects. Any systematic effects which may be due to position across the nar- 
row part of the beam will be neglected. 

CONSIDERATION OF RESIDUALS FROM THE MINITRACK NETWORK AS A WHOLE 

The Minitrack residuals were calculated on the basis of the 5-1/4-day standard orbit described 
under "Summary of Optical Data Used in the ’Standard Orbit’ Solution." In order to determine the 
best method of detecting possible systematic trends, the residuals were analyzed by considering 
the various possible methods of combining At and Am. Histograms for the entire set of residuals 
from all 12 stations are shown in Figure 6, with At and Am together and separately. There is no ob- 
vious evidence from these histograms to suggest any strong systematic differences between the op- 
tical data and the Minitrack data considered as a whole. No distinction has been made at this point 
according to the Minitrack mode of tracking (i.e., equatorial or polar). 

The residuals which were the basis of the histograms depicted in Figure 6 were then 
grouped according to the mode of tracking. These results are displayed in Figure 7 in the form of 


*See also Simas, V. R., and Watkins, E. R., Jr., as cited in earlier footnotes. 
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three-dimensional histograms showing relative frequency of magnitudes of the residuals in steps 
of 0.1 x 10“ 3 and by elevation angle in steps of 20°. The histogram labeled "Beam Only Effects" 
is obtained from the frequency of At eq andAm polar . The histogram labeled "Beam Only + Refraction 
Effects" is obtained from the frequency of At- . and Am residuals. The mean values and stand- 
ard deviations of the residuals in the elevation dimension are printed at the top. For the "Beam 
Only Effects," the mean values in the different elevation dimensions remain fairly close to 0. In 
the "Beam Only + Refraction Effects" histogram, there is a definite shift of the mean values of the 
residuals to the negative at low elevations. In the 10° to 30° elevation dimension, the shift amounts 
to -0.2 x 10 “ 3 . The small magnitude of the mean values of the "Beam Only Effects" set of resid- 
uals is further verification of the absence of strong systematic differences between optical and 
Minitrack data sets. The relatively large values of the standard deviations can be explained on the 
basis of individual Minitrack station biases. This effect will be shown in the next section. 

SYSTEMATIC ERRORS AT INDIVIDUAL MINITRACK SITES 

In order to detect possible biases at individual Minitrack station sites, the residuals were re- 
plotted on a station basis. Figures 8a through 81 depict the residuals for each station separately. 

In these plots by individual station, the l and m direction cosine residuals have been plotted sep- 
arately as a function of elevation angle. A circle represents a residual obtained from an observa- 
tion when the equipment was operating in the polar tracking mode. A dot represents a residual 
obtained from an observation when the equipment was operating in the equatorial mode. This same 
residual set has been summarized in histogram form in Figures 9a and 9b without regard to track- 
ing mode. There are pronounced biases, particularly at the stations of College, Alaska (COLEGE), 
and Mojave, California (MOJAVE). In particular, College has no positive residuals in At and no 
negative residuals in Am. At Mojave, there is strong positive bias in the At residuals and a strong 
negative bias in the Am residuals. 

Considerable caution must be exercised in attempting to analyze residual sets which may or 
may not contain biases. Also a certain ingenuity is often required in order to unmask and identify 
these biases if they do exist. For example, the histograms in Figure 6, which summarize the fre- 
quency of the Minitrack residuals from all stations both collectively (i.e., Art and Am together) and 
for At and Am separately, show no obvious biased tendency even in the separated Art and Am histo- 
grams. A rather broad dispersion is evident in this figure, however, both for the case where A i 
and Am are depicted together and where they are depicted separately. Using a different representa- 
tion (see Figure 9a), the residuals for some stations now show a strong biased tendency. In this 
case, the residuals have been separated and depicted by individual station as well as separated 
into the At and Am histograms. In this form of presentation, the biases at Mojave and College are 
pronounced. However, by considering the individual dispersions by At only and Am only for each 
station, the dispersions are in general found to be narrower than indicated by Figure 6 (histogram 
of all stations together). This suggests that if the sources of these systematic errors can be found 
and eliminated, there should be a corresponding reduction in the standard deviations both for the 
individual stations and for the collective data set from all stations. 
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ELEVATION VS DIRECTION COSINE RESIDUAL x 10 3 
b. Mojave, California 36 observations 
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Figure 8— Plots of elevation vs direction cosine residuals (elevation vs M x 10 3 
and elevation vs Am x 10 3 ) and mode of tracking (polar or equatorial). 
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Figure 8— Plots of elevation vs direction cosine residuals (elevation vs At x 10 3 
and elevation vs Am x 10 3 ) and mode of tracking (polar or equatorial). 
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Figure 8— Plots of elevation vs direction cosine residuals (elevation vs At x 10 3 
and elevation vs Am x 10 3 ) and mode of tracking (polar or equatorial). 
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Figure 8— Plots of elevation vs direction cosine residuals (elevation vs At x 10 3 
and elevation vs Am x 10 3 ) and mode of tracking (polar or equatorial). 
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Figure 8— Plots of elevation vs direction cosine residuals (elevation vs M x 10 3 
and elevation vs Am x 10 3 ) and mode of tracking (polar or equatorial). 
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Figure 8— Plots of elevation vs direction cosine residuals (elevation vs AC x 10 3 
and elevation vs Am x ]0 3 ) and mode of tracking (polar or equatorial). 
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RESIDUALS SEPARATED BY STATION AND INTO AND m 
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FTMYRS 


Figure 9a— Histogram of Minitrack residuals calculated from optically-determined 
"standard orbit"— residuals separated by station and into AfandAm. 


To illustrate the possibility of reducing the standard deviation of the data, a sample calcula- 
tion was made for Mojave (Figures 8b and 9a). Considering tJL and Am both separately and together, 
the following mean values and standard deviations were obtained: 



X 

M only 

+0. 18 x 

Am only 

-0.16 x 

Af + Am 

+0.01 x 


cr 


10~ 3 

0.18 x lO" 3 

10' 3 

0.05 x 10~ 3 

10- 3 

0.23 x 10’ 3 


It is seen that the standard deviation of At only is approximately 25 percent smaller than the stand- 
ard deviation of At and Am considered together. Next, the assumption was made that the biases for 
Mojave were known and that their effects caused the mean offsets of +0.18 x 10' 3 and -0.16 x 10~ 3 
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Figure 9b — Histogram of Minitrack residuals calculated from optically-determined 
’’standard orbit" — residuals separated by station and into M andAm. 


respectively in A t and Am. A further assumption was made that the effect of removing known biases 
could be approximated by subtracting the mean offsets from the corresponding residual sets. The 
standard deviation of the entire combined data set (l and m together) was then recalculated after the 
offsets had been subtracted from the original residual set. The new value of the standard devia- 
tion was 0.15 x 10“ 3 compared to the original value of 0.23 x 10“ 3 , a reduction of approximately 
33 percent. Although Mojave represented one of the most biased of the data sets used in this anal- 
ysis, this simple procedure reduced the standard deviation for Mojave to a value below the overall 
rms of fit of 0.19 x 10“ 3 for the entire data set. A further inspection of the Mojave residuals in 
Figure 8b indicates also that this bias is not primarily due to refraction- dependent residuals as 
defined earlier, since the Am polar residuals in Figure 8b show a strong negative bias although they 
are refraction-independent. Despite the indications of station biases noted in this section, however, 
it is felt that additional data sets should be examined before any definitive inferences of station 
biases can be properly made, since the data set analyzed extended for only 5-1/4 days. In a com- 
puter program which employs a rejection criterion for the use of data in orbital solutions, these 
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apparent station biases could cause an unwarranted rejection of data which might be basically good 
data. For example, it was noted in the fourth section (''Summary of Minitrack Data") that 10 per- 
cent of the data from College, Alaska, had residuals which exceeded 0.5 x 10~ 3 (the rejection 
criterion employed for analysis purposes in this paper). 


DIFFERENCES IN ORBITS ADJUSTED USING 
MINITRACK DATA ONLY AND OPTICAL DATA ONLY 

In order to obtain an estimate of the actual differences of position in inertial space as obtained 
by a Minitrack-determined orbit and an optically- determined orbit, trajectories obtained from the 
two adjusted solutions were calculated, differenced, and resolved into along- track, cross-track, and 
radial differences. The rms of fit of the orbit adjusted by Minitrack data only, not corrected for 
refraction or other effects, was 0.19 x 1CT 3 in terms of direction cosine for the 5-1/4-day arc used 
in the "standard" orbit solution. The position and velocity vectors for the epoch 01 H 38 M 22?000 
12/31/1965 UTC, which were adjusted on the basis of the 436 Minitrack observations, are given in 
Table 6. The Minitrack data included observations down to 15° elevation. The rms of fit of the 
standard orbit adjusted by optical data only was 3I T 0. In both this calculation and the previous anal- 
ysis, the start and end points of the Minitrack data were chosen to correspond to the equivalent 
start and end times of the optical data so that there would be no systematic orbital shift caused by 
overlap effects. The along-track, cross-track, and radial position differences of the first and last 
four hours of the Minitrack and optical orbits are shown in Figure 10. The along-track difference 
takes the approximate form of a sine curve with a two-hour period (the period of GEOS I) with an 
amplitude of approximately 110 meters superimposed upon a small secular term. The secular por- 
tion of the curve has a rate of approximately 16 meters/day. Considering the optical orbit as an 
error-free "reference orbit,” the root mean square of the total position errors over the 5-1/4-day 


Table 6 

Final Adjusted Position and Velocity Vectors at 
Epoch for the Minitrack-Determined Orbit. 


Minitrack-Adjusted 
Position Vector 


X: +5,690,533.8 meters 
Y: +1,474,657.9 meters 
Z: +6,013,372.3 meters 


Epoch: 

No. of Observations: 436 


M init r a ck-Ad justed 
Velocity Vector 

X: -4,685.6272 meters/sec. 
Y: +3,849.5020 meters/sec. 
Z: +2,939.1288 meters/sec. 


01 H 38 M 22SQ00 12/31/1965 UTC 


Arc Length: 5-1/4 days 

Rms of Fit: 1.906 x 10“ 4 
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ALONG -TRACK DIFFERENCE 

CROSS -TRACK DIFFERENCE 

x x x x RADIAL DIFFERENCE 



Figure 10— Differences between 5-1/4-day trajectories adjusted by Minitrack data only and optical data 
only using the SAO M -1 gravity model. (Differences are resolved along-track, cross-track and radially.) 


arc of the Minitrack orbit was approximately 165 meters. The Minitrack solution utilized data down 
to approximately 15° elevation, was not corrected for tropospheric or ionospheric refraction, and was 
not weighted according to elevation to compensate for any deterioration of the data at the lower 
elevations. 


CONCLUSIONS 

It has been shown that an optically-determined orbit is sufficiently accurate for evaluating 
Minitrack data. The standard deviation of position of the Minitrack-determined orbit analyzed in 
this report, using an optic ally- determined orbit as a standard, was approximately 165 meters for 
a 5-1/4-day arc. The Minitrack -determined orbit used data acquired down to 15° elevation. These 
data were not corrected for tropospheric or ionospheric refraction and were not weighted in the 
solution for elevation effects. The comparatively small standard deviation of 165 meters which 
resulted may be attributed in large part to a careful selection of computing tools used in the anal- 
ysis. In particular, the use of an extremely accurate gravity model (SAO M-l model) and a track- 
ing complement referred to a datum consistent with this gravity model (SAO C-5) contributed to 
the accuracy of the results obtained. (See Reference 16 for a more detailed analysis of gravity 
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effects.) The relatively small standard deviation also suggests that low-elevation Minitrack data 
when corrected for refraction effects may possibly be used with confidence in orbital solutions. 

The use of optical data as an independent data set revealed systematic offsets of the residuals 
at certain Minitrack sites. The use of such an independent data set could be effectively used to 
determine the calibration coefficients for the individual Minitrack stations and to investigate the 
feasibility of using the Minitrack data for self- calibration purposes. 

Goddard Space Flight Center 
National Aeronautics and Space Administration 
Greenbelt, Maryland, January 6, 1969 
311-07-21-01-51 
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Appendix A 

Brief Description of the 136-Mc Minitrack Interferometer Tracking System 


The basic Minitrack configuration consists of an antenna field of 13 antennas arranged as shown 
in Figure Al. In addition to the antenna system, there is an operations building housing the electronic 
systems. A MOTS 40-inch-focal -length camera is located at the exact center of the fine-beam ar- 
rays of the antenna system. This camera is normally used for calibrating the antenna system by 
photographing an airborne flashing light beacon against a background of stars while the antenna sys- 
tem is simultaneously receiving a 136-Mc radio signal from the same airborne source. Since theo- 
retically the physical location of the MOTS camera and the antenna system center are identical, a 
unique system is available for intercomparing the two distinct types of tracking systems. 
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The electronics of the Minitrack interferometer system is designed to locate satellites trans- 
mitting on 1000 separate frequencies ranging from 136.000 Me to 136.999 Me in 1-kc steps. The 
tracking is performed in a completely passive mode. The possibility of extending this range to 
138 Me is currently being considered. 

The antenna system is designed to be highly directive, resulting in a high signal gain. There 
are two possible modes of tracking, the T 'equatorial” mode and the "polar” mode. The wavelength 
dimensions displayed in Figure A1 are based upon a received signal whose frequency is 136.555 Me. 

The antenna arrays located at the ends of the north-south and east-west baselines are used 
in sets of four for each of the two tracking modes and are abbreviated as follows: 

- . Group 1 — Equatori al Tracking Mode 


Equatorial Fine-Beam Arrays 


NFE 

North Fine Equatorial 

SFE 

South Fine Equatorial 

EFE 

East Fine Equatorial 

WFE 

West Fine Equatorial 

Group 2 — 

Polar Tracking Mode 

Polar 

■ Fine- Beam Arrays 

NFP 

North Fine Polar 

SFP 

South Fine Polar 

EFP 

East Fine Polar 

WFP 

West Fine Polar 


The distance between NFE and SFE is the same as indicated between EFE and WFE. Also, the dis- 
tance between NFP and SFP is the same as indicated between EFP and WFP. 

The reception pattern of Group 1 (Equatorial Fine- Beam) is shown on the left side of Figure A2. 
The reception pattern of Group 2 (Polar Fine- Beam) is shown on the right side of the same figure. 
These two groups comprise the accurate portion of the antenna system. The individual stations can 
electronically switch between the polar and equatorial modes of tracking very quickly, so that both 
modes of tracking are possible on the same pass of a satellite. The actual shape of the wedge or 
fan-shaped pattern of the fine beam is defined by the locus of points where the signal from the satel- 
lite's Minitrack beacon drops 3 decibels from the zenith signal. Under normal conditions, the sat- 
ellite is tracked only when it is located with this wedge-shaped fine beam. Tracking is usually 
further restricted to a maximum distance of 20° each way from the zenith in the long-beam direction 
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ZENITH 



The Minitrack Observations are 
1 ( = cos a) 
and m ( = cos/9) 


ZENITH 


EQUATORIAL TRACKING MODE 
( Long beam oriented N - S ) 

A t - Dependent upon long beam only 

Am - Dependent upon long beam & refraction 


"v.E 



POLAR TRACKING MODE 
( Long beam oriented E -W) 

A t - Dependent upon long beam & refraction 
Am - Dependent upon long beam only 


Figure A2 — Approximate reception pattern of the fine beam of the 136-Mc Minitrack antenna array. 


of this pattern in order to minimize refraction 
effects, A typical reception pattern for the 
EFE array in the east- west direction at the 
Minitrack station at Santiago, Chile, is shown 
in Figure A3. 

In order to calculate the direction cosines 
of a satellite, it is necessary to know the ab- 
solute phase difference of the signal as re- 
ceived simultaneously by a pair of the fine- 
beam antennas, for example, WFE and EFE. 

Since only relative phase differences of the 
signal can be measured as received by the 
fine-beam pairs, it is necessary to introduce 
five additional arrays to resolve these phase-difference ambiguities which arise. The ambiguity 
antennas are labeled in Figure A1 as W, N, C, E, S (West, North, Common, East, and South respec- 
tively). This group of five antennas can be used to resolve the ambiguities for both the Group 1 
and Group 2 fine-beam arrays. The W-N and E-S combinations are used for east-west medium and 
north-south medium resolution respectively. The N-C and C-E combinations are used for north- 
south and east-west coarse resolution respectively. In being so used, the combinations are sep- 
arated by 4 and 3.5 wavelengths respectively, so that the difference in separation of the medium 


NOTE: THIS FIGURE IS REPRODUCED FROM 
"INSPECTION OF THE 136.5 MC 
MINITRACK STADAN ANTENNAS AT 
SANTIAGO, CHILE, S.A.," 
JANUARY 1965, NEW MEXICO 
STATE UNIV. 



30 40 50 

DEGREES FROM CENTER 


Figure A3— Antenna pattern for east fine 
equatorial Minitrack antenna at Santiago, Chile. 
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and coarse combinations used to resolve either north-south or east-west ambiguities is only 1/2 
of a wavelength. 


SATELLITE 


The electronics of the Minitrack system divides the 360 electrical degrees phase measurement 
into 1000 parts of 0.36 electrical degrees each. This measurement represents the upper precision 

limit of the equipment. If it is assumed that 
the fine-beam baseline is approximately 100 
meters (see Figure A4), the resulting accuracy 
in terms of space angle can be estimated ap- 
proximately by the following procedure. The 
direction cosine l (considering the east- west 
baseline) is related to the electrical phase 
angle by the equation 



cos a - o 


n\ 


_ B = 100 Meters_ 
( Baseline ) 


Figure A4— Simp! ified geometry of satellite 136-Mc 
Minitrack beacon signal as received by the west fine 
equatorial and east fine equatorial antenna arrays. 

If 0 is measured in degrees, the equation becomes 

't = cos a - 

To find the error in l and a, 


_0 k 

277 B 

where 

0 = electrical angle in radians, 

A = signal wavelength, 

= 2.2 meters at 136.555 Me, 

B = antenna separation distance or base- 
line length, 

= 100 meters, 

d = phase difference in linear measure, 

n = phase difference in number of wave- 
lengths of signal. 


0 2 . 2 
360° 100 


zve 


- sin aAa 


2.2 

360 x 100 


A0 


2 . 2 
3.6 


X 10" 4 M> • 


32 



If = 0.36 degrees, then 


- - sinaAa - ^ * 0.36 x 10 4 

= 2 . 2 x 10" 5 • 

For a near-overhead pass, sin a ~ 1. Neglecting the minus sign, 

Act = 2. 2 x 10~ 5 radians 

= 2. 2 x 10“ 5 x (2 x 10 5 )" 

= 4'M - 
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Appendix B 


Preprocessing of Optical Observations 
Preprocessing of Optical Data 

The first step in the processing of optical observations (and one usually performed by the ob- 
serving source) consists of developing a plate or film and identifying thereon the image or images 
of the satellite and the images of several reference stars whose right ascensions and declinations 
are well known. The initial measurements both of satellite images and of reference stars consist 
of linear rectangular coordinates. From the knowledge of the spherical coordinates of the refer- 
ence stars, the right ascensions and declinations of the satellite images may be calculated. These 
coordinates as received by the preprocessor may be referred to the mean equator and equinox of 
date, true equator and equinox of date, or mean equator and equinox of some standard epoch. 

The preprocessor then transforms these observations to a common coordinate system. It 
transforms all right ascensions and declinations to the true equator and equinox of the epoch of the 
observations being processed. If the observations were originally referred to the mean equator 
and equinox of a particular epoch, it is necessary only to precess from that epoch to the epoch of 
the observations. However, if they were referred to the true equator and equinox of a particular 
epoch, it is necessary first to transform them to the mean equator and equinox of that same epoch 
and then precess to the epoch of the observations. Finally, a transformation must be made from 
the mean equator and equinox of the epoch of the observations to the true equator and equinox of 
the epoch of the observations. 

Nutation 

The transformation from the true equator and equinox of date to the mean equator and equinox 
of date is 


Y = NX , 

where 



cos 

8 

cos a 



m 

m 

= 

cos 

8 

sin a 



m 

m 


sin 





cos 

S T 

cos a T 

= 

cos 

S T 

sin a T 


sin 

S T 
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and 


N 


A0 COS € 

' m 

1 

“Ae 


A^A sin e 

T m 

Ae 


“A ip cos e a 
“Av// sine 


1 


a m , S m = right ascension and declination referred to mean equator and equinox of date, 
a T , S T = right ascension and declination referred to true equator and equinox of date, 
e m = mean obliquity of date, 

A p = nutation in longitude, and 
Ae = nutation in obliquity. 

The inverse transformation is simply: 


X = N ~ 1 Y = N t Y . 


Precession 

The transformation from the mean equator and equinox of 1950.0 to the mean equator and 
equinox of an arbitrary epoch ti is 


where 


Y = px , 



cos S ti 

cos a fcl 


= 

C °s 8 tl 

sina «! 

■ 


_ sin 5 ti 

__ 


cos 

s 

1950-0 

cos a 1 

1950-0 

cos 

^ 1950-0 

Sina i9S0 

.0 

sin 

8 

1950-0 




and 


P 


( cos z cos 0 cos £ - sin z sin £ ) 
( sin z cos 0 cos £ + cos z sin £) 
( sin 6 cos £) 


(- cos z cos 6 sin £ - sin z cos £) (- cos z sin 0) 

(“ sin z cos 0 sin £ + cos z cos £) (- sin z sin 6 ) 

(-sin0sin£) (cos 0) 
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The inverse transformation is 


X = P ' 1 Y = P T Y . 

Since the expressions in z, i are tied to 1950.0 as an epoch, the precession between two 
different epochs, neither of which is 1950.0, must be performed in two steps, using 1950.0 as an 
intermediary epoch. The above expression for p is rigorous. There are simple third-degree 
polynomials in time derived by Newcomb which permit the calculation of z, 0, and i. There exists 
an even simpler form of the matrix p which permits the calculation of its elements by means of 
third-degree polynomials expressed directly in terms of the variable t (time). This simplification 
bypasses the necessity of calculating the sines and cosines of the angles z, e 9 and £. These sim- 
plified matrix elements are derived by expanding the sines and cosines of z, 0 9 and £, contained in 
the elements of p, into a series, performing the necessary multiplications, and dropping terms 
exceeding the third degree. The appropriate polynomial expressions in t are then substituted into 
the remaining expressions containing z, <9, and 5. After the necessary multiplications have again 
been performed, all terms in t higher than the third degree are dropped. The final expression 
for p then consists simply of nine elements in terms of a third-degree polynomial in time. 
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Appendix C 

Force Models Used in NONAME 


Force Models 

The data reduction program in its present form incorporates four force models. These are: 

1. The earth's gravitational field; 

2. The solar and lunar gravitational perturbations; 

3. Solar radiation pressure; 

4 . Atmo spheric drag. 

The program is so designed that the gravitational coefficients and pertinent physical character- 
istics of satellites, such as reflectivity, cross-sectional area mass, and drag coefficient, can be 
simply changed through card input or block data statement. 


The Earth’s Gravitational Field 

The formulation of the geopotential used is: 


u 



EE®' 


n = 2 m= 0 


P m ( sin <p) \ C cos mk + S 

n ' ^ ' L nm nm 



(Cl) 


where 


G 

M 

r 

a 

k 

P n m ( sin 0) 


is the universal gravitational constant, 

is the mass of the earth, 

is the geocentric distance to the satellite 

is the earth's mean equatorial radius, 

is the sub- satellite latitude, 

is the sub- satellite east longitude, and 

is the associated spherical harmonic of degree n and order m. 


The design of the potential function requires that normalized gravitational coefficients c n m and 
S n m be used. The program is presently capable of accepting coefficients up to (20, 20) or any sub- 
set of these. 


The SAO M-l earth gravitational model (Reference Cl) modified by the GEOS I resonant har- 

13, 12* ^ 14 , 12 * ‘ J 14,12 j ^15. 12* 


monies (c„ s,, C,„ C 1S S 15 ,,) (Reference C2) is listed in Table Cl. These 
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Table Cl 


SAO M-l Harmonic Coefficients (Normalized) . 














n 

m 

C X 10 6 

S X 10 6 

n 

m 

C xio 6 

S x 10 6 

2 

0 

-484.1735 


8 

2 

0.026 

0.039 

2 

1 



8 

3 

-0.037 

0.004 

2 

2 

2.379 

-1.351 

8 

4 

-0.212 

-0.012 





8 

5 

-0.053 

0.118 

3 

0 

0.9623 


8 

6 

-0.017 

0.318 

3 

1 

1.936 

0.266 

8 

7 

-0.0087 

0.031 

3 

2 

0.734 

-0.538 

8 

8 

-0.248 

0.102 

3 

3 

0.561 

1.620 

9 

0 

0.0122 


4 

0 

0.5497 


9 

1 

0.117 

0.012 

4 

1 

-0.572 

-0.469 

9 

2 

-0.0040 

0.035 

4 

2 

0.330 

0.661 





4 

3 

0.851 

-0.190 

10 

00 

0.0118 


4 

4 

-0.053 

0.230 

10 

01 

0.105 

-0.126 





10 

02 

-0.105 

-0.042 

5 

0 

0.0633 


10 

03 

-0.065 

0.030 

5 

1 

-0.079 

-0.103 

10 

04 

-0.074 

-0.111 

5 

2 

0.631 

-0.232 





5 

3 

-0.520 

0.007 

11 

00 

-0.0630 


5 

4 

-0.265 

0.064 

11 

01 

-0.053 

0.015 

5 

5 

0.156 

-0.592 

12 

00 

0.0714 


6 

0 

-0.1792 


12 

01 

-0.163 

-0.071 

6 

1 

-0.047 

-0.027 

12 

02 

-0.103 

-0.0051 

6 

2 

0.069 

-0.366 

12 

12 

-0.031 

0,0008 

6 

3 

-0.054 

0.031 

13 

00 

0.0219 


6 

4 

-0.044 

-0.518 

13 

12 

-0.06769 

0.06245 

6 

5 

-0.313 

-0.458 

13 

13 

-0.059 

0.077 

6 

6 

-0.040 

-0.155 





7 

0 

0.0860 . 


14 

00 

-0.0332 


7 

1 

0.197 

0.156 

14 

01 

-0.015 

0.0053 

7 

2 

0.364 

0.163 

14 

11 

0.0002 

-0.0001 

7 

3 

0.250 

0.018 

14 

12 

0.00261 

-0.02457 

7 

4 

-0.152 

-0.102 

14 

14 

-0.014 

-0.003 

7 

5 

0.076 

0.054 





7 

6 

-0.209 

0.063 

15 

09 

-0.0009 

-0.0018 

7 

7 

0.055 

0.096 

15 

12 

-0.07473 

-0.01026 

8 

0 

0.0655 


15 

13 

-0.058 

-0.046 

8 

1 

-0.075 

0.065 

15 

14* 

0.0043 

-0.0211 


The normalized coefficients (c , s ) are related to the denormalized coefficients (c , s ) as indicated 

' n.m* n.nu 'n.m'n.ra' 

below: 


C(m, m) = £(n-m)!(2n + l )k/(n + m) ! ] 1//2 C(n , m) 
S(n, m) = £(n - m) ! (2n + l)k/(n + m) ! ] 1/2 S(n, m) 

k - l when m = o 
k = 2 when m ¥ o 
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coefficients have been used extensively in the NONAME orbit determination program for the re- 
duction of GEOS I optical and electronic data. The same data sets have been reduced using various 
other gravity models. An inter comparison of the results can be found in Reference C3. 

The transformation of the geopotential in earth- fixed coordinates (r, <f> 9 k) to gravitational ac- 
celerations in inertial coordinates (x, y, z) is accomplished as follows: 

du du d_/V . 

x © " dr dx + Tcf) dx + dk dx (02) 


and similarly for y 0 , z 0 , where the subscript denotes accelerations caused by the earth’s field. 


Solar and Lunar Gravitational Perturbations 


The perturbations caused by a third body, e.g., the sun or moon, on a satellite orbit are treated 
by defining a disturbing function r (Reference C4) which can be treated as the potential function u. 
The solar perturbation R e takes the form 


where 


GMm^ 


/ 2r r 2 V 1/2 rs 

-n 


(C3) 


s = 


( ? - 7 o) 


and 

m Q is the mass of the sun in earth masses, 
r 0 is the geocentric position vector of the sun, 
r Q is the geocentric distance to the sun, 
r is the geocentric position vector of the satellite, 
r is the geocentric distance to the satellite, 

G is the universal gravitational constant, and 
M is the mass of the earth. 

The satellite acceleration due to the sun is then 

d R 0 d_r dtf> dk . 

x © " dr dx + d<P dx + dk dx (C4) 
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■plus y Q , z 0 , similarly derived, where A. and <f> are the longitude and latitude of the satellite respec- 
tively. The lunar perturbation is found from Equation C3 by substituting the lunar mass and dis- 
tance for those of the sun. 

The lunar and solar ephemerides are computed internal to the program. These positions are 
computed at 10 equal intervals over each five-day period and least squares fit to a fourth-order 
polynomial in time about the midpoint of the five-day period. The positions of these bodies are then 
determined at each data point by evaluating the polynomial at the observation time. 

Solar Radiation Pressure 

The satellite acceleration due to solar radiation pressure is formulated as follows (Refer- 
ence C5): 


RAD 


J v L ; 

m ' x * 


(C5) 


plus y RAD , z RAD , similarly derived, where 

L is the inertial unit vector from the geocenter to the sun and has the components L x , L y , 

L z j 

A is the cross-sectional area of the satellite, 
m is the satellite mass, 

y is a factor depending on the reflective characteristics of the satellite, 

v is the eclipse factor such that: 

fO when satellite is in earth's shadow 

v - A 

U when satellite is illuminated by the sun, and 
P Q is the solar radiation pressure in the vicinity of the earth, 


4.5 X 10" 6 


Newton 


At present, it is assumed that the satellite is specularly reflecting with reflectivity p 9 and thus 

r = (i +p) • (c 6 ) 

The vector L and the eclipse factor z> are determined from the solar ephemeris subroutine 
previously described, and from the satellite ephemeris, and are based upon the approximation of 
a cylindrical earth shadow. 
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Atmospheric Drag 

The atmospheric decelerations are computed as follows: 

^D AVV x . 

X DRAG 2m ’ (C7) 


plus y DRAG , z DRAG , similarly derived, where 
P is the ambient atmospheric density, 

C D is the satellite drag coefficient, 

A is the projected area of the satellite on a plane perpendicular to direction of motion, 
m is the satellite mass. 

The velocity vector v, given in inertial coordinates by 

v - vi+vj+vlc, 

x y J z (C8) 

can be chosen to be either the velocity relative to the atmosphere, which implies that the atmosphere 
rotates with the earth, or the inertial velocity, which assumes that the atmosphere is static. Pres- 
ently, the former assumption is made. 

The density p is computed from the 1962 U. S. Standard Atmosphere. 
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Appendix D 


Tracking Station Coordinates 
Datum Parameters and Station Coordinates 

For the purpose of long-arc satellite data reduction and intercomparison, all GEOS I partici- 
pating tracking stations have been transformed to a common datum. The common datum selected 
is the SAO Standard Earth C-5 model (semi-major axis = 6378165 meters, flattening coefficient 
= 298.25) (Reference Dl)in which the Baker -Nunn stations are used as the controlling stations for 
all other stations to be transformed. Descriptions and formulations to effect the transformations 
from major and isolated datums are presented in Reference D2. The transformation of local datum 
station coordinates to a common "center-of-mass reference" system is important to accomplish, 
since the datum shifts are quite large. For example, on the North American Datum the center-of-mass 
shift to the C-5 Standard Earth is approximately 250 meters. The center-of-mass coordinates of the 
SAO C-5 Baker-Nunn stations are assessed by SAO to have approximately 20-meter accuracy. 

In order to effect any transformation, the parameters of the original datums must be known as 
well as the geodetic latitude, longitude, and height. Table D1 provides a listing of the original 
datums and their parameters on which the stations were originally surveyed. Tables D2 through 
Dll list alternately the original surveyed ellipsoidal positions and the SAO C-5 ellipsoidal positions 
for over 100 GEOS I tracking stations that have been used in the long-arc inter comparison effort. 
These tables contain symbols designating the source of original coordinates. The symbols are de- 
fined in the following section, with a list of sources. The C-5 positions for 1TANAN and MADGAR 
(Reference D3) have been derived by the station estimation technique contained in the Orbit Determi- 
nation Program NONAME. Tables D12 through D21 provide a listing of the station names from 
which the six- letter designations have been derived. 

Sources 

The following sources were used to obtain the original datum positions: 

Symbol Source 

A Geodetic Parameters for a 1966 Smithsonian Institu- 

tion Standard Earth (Reference Dl). 

B Goddard Directory of Tracking Station Locations; 

August 1966; Goddard Space Flight Center. 

C NWL-8 Geodetic Parameters Based on Doppler 

Satellite Observations; July 1967; R. Anderle and 
S. Smith, Naval Weapons Laboratory. 
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Since the above official documents did not contain all those positions that were to be trans- 
formed, it was necessary to contact other sources for the positions of the remaining stations. 
These sources are: 

Symbol Source 

D Private communication with personnel at SAO 

(K. Haramundanis, B. Miller, A. Girnius). 

E Private communication with 1381st Geodetic Survey 

Squadron, USAF (S. Tischler). 

F Private communication with personnel at USC&GS 

(B. Stevens). 

G Private communication with personnel at U. S. Army 

Engineers Topographic Laboratories (L. Gambino). 

H Private communication with NASA Space Science 

Data Center (J. Johns, D. Tidwell). 

I General Station Data Sheet— GEOS A Project Manager, 

NASA Headquarters. 


Table D1 


Parameters of Original Datums. 


Datum Name 

Semi-Major Axis 
(meters) 

i/f 

North American (NA.) 

6378206.4 

294.9787 

European 

6378388.0 

297.0 

Tokyo 

6377397.2 

299.1528 

Argentina 

6378388.0 

297.0 

Mercury 

6378166.0 

298.3 

Madagascar 

6378388.0 

297.0 

Australian Nat'l. 

6378160.0 

298.25 

Old Hawaiian 

6378206.4 

294.9787 

Indian 

6377276.3 

300.8017 

Arc (Cape) 

6378249.1 

293.4663 

1966 Canton ASTRO 

6378388.0 

297.0 

Johnston Island 1961 

6378388.0 

297.0 

Midway ASTRO 1961 

6378388.0 

297.0 

Navy IBEN ASTRO 1947 

6378206.4 

294.9787 

Provisional DOS 
ASTRO 1962, 65 

6378388.0 

297.0 

Allen Sodano Lt, 

6378388.0 

297.0 

1966 SECOR ASTRO 

6378388.0 

297.0 

Viti Levu 1916 

6378249.1 

293.4663 

Corrego Alegre 

6378206.4 

294.9787 

USGS 1962 ASTRO 

6378206.4 

294.9787 

Berne 

6377397.2 

299.1528 
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Table D2 


SAO— Optical— Source A. 







Geodetic 


Source 

Name 

Station No. 

Latitude 

Longitude 

Height 

(meters) 

Datum 


1 ORGAN 

9001 

32°25'24V56 

253°26'51'.'17 

1649 

N.A. 




32 25 24.70 

253 26 48.29 

1610 

C-5 


lOLFAN 

9002 

-25 57 33.85 

28 14 53.91 

1562 

Arc (Cape) 




-25 57 37.67 

28 14 51.45 

1560 

C-5 


WOOMER 

9003 

-31 06 07.26 

136 46 58.70 

162 

Australian 




-31 06 04.14 

136 47 01.93 

158 

C-5 


IS PAIN 

9004 

36 27 51.24 

353 47 41.47 

7 

European 




36 27 46.68 

353 47 36.55 

56 

C-5 


1 TOKYO 

9005 

35 40 11.08 

139 32 28.22 

58 

Tokyo 




35 40 23.03 

139 32 16.42 

84 

C-5 


1NATOL 

9006 

29 21 38.90 

79 27 25.61 

1847 

European 




29 21 34.38 

79 27 27.05 

1855 

C-5 


1QUIPA 

9007 

-16 28 05.09 

288 30 22.84 

2600 

N.A. 




-16 27 58.04 

288 30 24.02 

2479 

C-5 


1SHRAZ 

9008 

29 38 17.90 

52 31 11.80 

1578 

European 




29 38 13.59 

52 31 11.20 

1561 

C-5 


1CURAC 

9009 

12 05 21.55 

291 09 42.55 

23 

N.A. 




12 05 24.93 

291 09 43.97 

-33 

C-5 


1JUPTR 

9010 

27 01 13.00 

279 53 12.92 

26 

N.A. 




27 01 14.23 

279 53 12.95 

-36 

C-5 


1VILDO 

9011 

-31 56 36.53 

294 53 39.82 

598 

Argentinean 




-31 56 36.35 

294 53 36.11 

636 

C-5 


1MAUIO 

9012 

20 42 37.49 

203 44 24.11 

3027 

Old Hawaiian 




20 42 25.66 

203 44 33.23 

3027 

C-5 


AUSBAK 

9023 

-31 23 30.82 

136 52 39.02 

141 

Australian 




-31 23 27.69 

136 52 42.23 

137 

C-5 


OSLONR 

9426 

60 12 40.38 

10 45 08.74 

585 

European 




60 12 38.88 

10 45 02.26 

573 

C-5 

I 

NATAL B* 

9029 

-05 55 50.00 

324 50 18.00 

112 

NA. 




-05 55 43.49 

324 50 21.30 

45 

C-5 

D 

AGASSI* 

9050 

42 30 20.97 

288 26 28.71 

193 

N.A. 




42 30 20.51 

288 26 29.79 

138 

C-5 

I 

COLDLK* 

9424 

54 44 38.02 

249 57 25.85 

597 

NA. 




54 44 37.26 

249 57 21.90 

548 

C-5 

I 

EDWAFB* 

9425 

34 57 50.68 

242 05 11.39 

784 

N.A. 




34 57 50.17 

242 05 07.80 

754 

C-5 

I 

RIGLAT* 

9428 

56 56 54.00 

24 03 42.00 

5 

European 




56 56 52.37 

24 03 37.49 

-15 

C-5 

I 

POTDAM* 

9429 

52 22 55.00 

13 04 01.00 

111 

European 




52 22 52.33 

13 03 55.80 

106 

C-5 

I 

ZVENIG* 

9430 

55 41 37.70 

36 46 03.00 

145 

European 




55 41 36.17 

36 46 00.17 

114 

C-5 


*These SAO station positions were derived by using the weighting scheme described in Reference D2 (in its second section, "Coordinate 
Transformation”). 
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Table D3 

STADAN — Optical — Source B.* 


Name 

Station No. 

Latitude 

Longitude 

Geodetic 

Height 

(meters) 

Datum 

1BPOIN 

1021 

38°25'49763 

282°54 T 48 r . r 23 

5 

NA. 



38 25 49.44 

282 54 48.65 

-50 

C-5 

1FTMYR 

1022 

26 32 51.89 

278 08 03.93 

19 

NA. 



26 32 53.08 

278 08 03.80 

-42 

C-5 

lOOMER 

1024 

-31 23 30.07 

136 52 11.05 

152 

Australian 



-31 23 26.96 

136 52 14.25 

148 

C-5 

1QUITO 

1025 

-0 37 28.00 

281 25 14.81 

3649 

NA. 



-0 37 22.63 

281 25 15.23 

3554 

C-5 

1LIMAP 

1026 

-11 46 44.43 

282 50 58.23 

155 

N A. 



-11 46 37.56 

282 50 58.86 

34 

C-5 

1SATAG 

1028 

-33 09 07.66 

289 19 51.35 

922 

N.A. 



-33 08 58.76 

289 19 52.59 

705 

C-5 

1MOJAV 

1030 

35 19 48.09 

243 06 02.73 

905 

NA. 



35 19 47.57 

243 05 59.18 

874 

C-5 

1JOBUR 

1031 

-25 52 58.86 

27 42 27.93 

1530 

Arc (Cape) 



-25 53 02.70 

27 42 25.41 

1546 

C-5 

1NEWFL 

1032 

47 44 29.74 

307 16 43.37 

104 

NA. 



47 44 28.73 

307 16 46.67 

58 

C-5 

1COLEG 

1033 

64 52 19.72 

212 09 47.17 

162 

N.A. 



64 52 17.78 

212 09 37.29 

139 

C-5 

1GFORK 

1034 

48 01 21.40 

262 59 21.56 

253 

NA. 



48 01 20.81 

262 59 19.55 

200 

C-5 

1WNKFL 

1035 

51 26 44.12 

359 18 14.62 

62 

European 



51 26 40.67 

359 18 08.35 

76 

C-5 

1ROSMA 

1042 

35 12 06.93 

277 07 41.01 

914 

NA. 



35 12 07.03 

277 07 40.81 

857 

C-5 

1TANAN 

1043 

-19 00 27.09 

47 18 00.46 

1377 

Tananarive 



-19 00 33.26 

47 17 58.89 

1355 

C-5 


*The coordinates in the above table are Identical to the corresponding Minitrack site coordinates if the latter exist, since the STADAN 
MOTS cameras are colocated with the center of the Minitrack antenna array for calibration purposes. 


Table D4 

STADAN—R/R— Source B. 






Geodetic 


Name 

Station No. 

Latitude 

Longitude 

Height 

Datum 





(meters) 


MADGAR 

1122 

-19 o 01' 13732 

47°18’09V45 

1403 

Tananarive 



-19 01 19.41 

47 18 07.96 

1382 

C-5 

ROSRAN 

1126 

35 11 45.05 

277 07 26.23 

880 

NA. 



35 11 45.15 

277 07 26.02 

823 

C-5 

CARVON 

1152 

-24 54 14.85 

113 42 55.05 

38 

Australian 



-24 54 12.29 

113 42 58.54 

10 

C-5 
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Table D5 


Navy TRANET— Doppler— Source C. 






Geodetic 


Name 

Station No. 

Latitude 

Longitude 

Height 

(meters) 

Datum 

LASHAM 

2006 

51°11'10V62 

358°58'30751 

182 

European 



51 11 07.12 

358 58 24.25 

196 

C-5 

SANHES 

2008 

-23 13 01.74 

314 07 50.59 

608* 

Corrego 

Alegre 



-23 13 01.74 

314 07 50.59 

608 

C-5 

PHILIP 

2011 

14 58 57.79 

120 04 25.98 

8 

Tokyo 



14 59 16.42 

120 04 21.61 

-70 

C-5 

SMTHFD 

2012 

-34 40 31.31 

138 39 12.39 

39 

Australian 



-34 40 28.16 

138 39 15.66 

31 

C-5 

MISAWA 

2013 

40 43 04.63 

141 20 04.69 

-10 

Tokyo 



40 43 14.63 

141 19 51.45 

38 

C-5 

ANCHOR 

2014 

61 17 01.98 

210 10 37.46 

61 

NA. 



61 16 59.60 

210 10 28.60 

44 

C-5 

TAFUNA 

2017 

-14 19 50.19 

189 17 13.96 

6* 

USGS 

1962 ASTRO 



-14 19 50.19 

189 17 13.96 

6 

C-5 

THULEG 

2018 

76 32 18.62 

291 13 46.72 

43 

N.A. 



76 32 20.72 

291 13 51.07 

-7 

C-5 

MCMRDO 

2019 

-77 50 51.00 

166 40 25.00 

-43 

Mercury 



-77 50 50.58 

166 40 35.02 

-29 

C-5 

WAHIWA 

2100 

21 31 26.86 

202 00 00.63 

380 

Old Hawaiian 



21 31 14.95 

202 00 09.83 

368 

C-5 

LACRES 

2103 

32 16 43.75 

253 14 48.25 

1201 

N.A. 



32 16 43.91 

253 14 45.34 

1162 

C-5 

LASHM2 

2106 

51 11 12.32 

358 58 30.21 

187 

European 



51 11 08.82 

358 58 23.95 

201 

C-5 

APLMND 

2111 

39 09 47.83 

283 06 11.07 

146 

NjV. 



39 09 47.59 

! 283 06 11.52 

90 

C-5 

PRETOR 

2115 

-25 56 46.09 

28 20 53.00 

1417 

European 



-25 56 49.97 

28 20 50.67 

1595 

C-5 

SHEMYA 

2739 

52 43 01.52 

174 06 51.43 

44 

NA. 



52 42 56.52 

174 06 44.17 

89 

C-5 

BELTSV 

2742 

39 01 39.46 

283 10 27.25 

50 

NJV. 



39 01 39.23 

283 10 27.72 

-5 

C-5 

STNVIL 

2745 

33 25 31.57 

269 09 10.70 

44 

NA. 



33 25 31.76 

269 09 09.66 

-10 

C-5 


*MSL. 
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Table D6 

Air Force— Optical — Source I.* 







Geodetic 


Source 

Name 

Station No. 

Latitude 

Longitude 

Height 

Datum 






(meters) 



ANTIGA 

3106 

17°08 T 5l”68 

298°12'37V41 

7 

NA. 




17 08 53.88 

298 12 39.19 

-42 

C-5 

E 

GRNVLE 

3333 

33 28 48.97 

268 59 49.17 

45 

N.A. 




33 28 49.15 

268 59 48.12 

-9 

C-5 


GRVILL 

3334 

33 25 31.95 

269 05 11.35 

43 

N.A. 




33 25 32.14 

269 05 10.30 

-10 

C-5 


USAFAC 

3400 

39 00 22.44 

255 07 01.01 

2191 

N.A. 




39 00 21.99 

255 06 58.32 

2147 

C-5 

E 

BEDFRD 

3401 

42 27 17.53 

288 43 35.03 

88 

N.A. 




42 27 17.06 

288 43 36.14 

33 

C-5 

E 

SEMMES 

3402 

30 46 49.35 

271 44 52.37 

79 

N.A. 




30 46 49.85 

271 44 51.64 

23 

C-5 


SWANIS 

3404 

17 24 16.57 

276 03 29.87 

83 

N.A. 




17 24 18.90 

276 03 29.71 

18 

C-5 


GRDTRK 

3405 

21 25 47.05 

288 51 14.03 

7 

NA. 




21 25 48.69 

288 51 15.03 

-48 

C-5 


CURACO 

3406 

12 05 22.11 

291 09 43.76 

23 

N.A. 




12 05 25.49 

291 09 45.16 

-34 

C-5 


TRNDAD 

3407 

10 44 32.78 

298 23 23.67 

269 

N.A. 




10 44 36.16 

298 23 25.43 

210 

C-5 


GRANFK 

3451 

47 56 38.63 

262 37 11.21 

296 

N.A. 




47 56 38.03 

262 37 09.15 

242 

C-5 


j TWINOK 

3452 

36 07 25.69 

262 47 04.48 

312 

N.A. 


I 


36 07 25.58 

262 47 02.68 

262 

C-5 


ROTHGR 

3453 

51 25 00.00 

9 30 06.00 

351 

European 




51 24 57.05 

9 30 00.58 

352 

C-5 


ATHNGR 

3463 

37 53 30.00 

23 44 30.00 

16 

European 




37 53 26.07 

23 44 26.73 

23 

C-5 


TORRSP 

3464 

40 29 18.53 

356 34 41.24 

588 

European 




40 29 14.10 

356 34 36.06 

635 

C-5 


CHOFUJ 

3465 

35 39 57.00 

139 32 12.00 

49 

Tokyo 




35 40 08.96 

139 32 00.19 

75 

C-5 


KINDLY 

3471 

32 22 57.30 

295 19 00.46 

26 

N .A. 




32 22 57.41 

295 19 02.09 

-23 

C-5 

E 

HUNTER 

3648 

32 00 05.87 

278 50 46.36 

17 

N.A. 




32 00 06.32 

278 50 46.32 

-40 

C-5 


JUPRAF 

3649 

27 01 14.80 

279 53 13.72 

26 

N.A. 




27 01 16.02 

279 53 13.72 

-37 

C-5 

E 

ABERDN 

3657 

39 28 18.97 

283 55 44.56 

4 

N.A. 




39 28 18.71 

283 55 45.10 

-51 

C-5 

E 

HOMEST 

3861 

25 30 24.69 

279 36 42.69 

18 

N.A. 




25 30 26.02 

279 36 42.70 

-44 

C-5 


CHYWYN 

3902 

41 07 59.20 

255 08 02.65 

1890 

N.A. 




41 07 58.61 

255_07 59.94 

1845 

C-5 


•Unless "Source” is specified otherwise. 
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Table D7 


Army Map Service— SECOR— Source H.* 


Source 

G 

I 

I 

I 

G 


Name 

HERNDN 

CUBCAL 

LARSON 

WRGTON 

GREENV 

TRUKIS 

SWALLO 

KUSAIE 


Station No. 

5001 

5200 

5201 

5202 
5333 

5401 

5402 

5403 


GIZZOO 

TARAWA 

NANDIS 

CANTON 

JONSTN 


5404 

5405 

5406 

5407 

5408 


MIDWAY 


5410 


G 

G 

G 


MAUIHI 

FTWART 

HNTAFB 

HOMEFL 


5411 

5648 

5649 
5861 


*Unless "Source” is specified otherwise. 
'MSL. 




Geodetic 


Latitude 

Longitude 

Height 

(meters) 

Datum 

38°59'37V69 

282°40 ! 16768 

119 

NA. 

38 59 37.47 

282 40 17.08 

64 

C-5 

32 48 00.00 

242 52 00.00 

101 

NA. 

32 47 59.74 

242 51 56.55 

71 

C-5 

47 11 00.00 

240 40 00.00 

354 

NA. 

47 10 58.76 

240 39 55.68 

319 

C-5 

43 39 00.00 

264 25 00.00 

481 

NA. 

43 38 59.49 

264 24 58.27 

428 

C-5 

33 25 32.34 

269 05 10.78 

43 

NA. 

33 25 32.53 

269 05 09.73 

-10 

C-5 

7 27 39.30 

151 50 31.28 

5* 

Navy IBEN 
ASTRO 1947 

7 27 39.30 

151 50 31.28 

5 

C-5 

-10 18 21.42 

166 17 56.79 

9 f 

1966 SECOR 
ASTRO 

-10 18 21.42 

166 17 56.79 

9 

C-5 

5 17 44.43 

163 01 29.88 

7 1 " 

ASTRO 1962, 
65 Allen 
Sodano Lt. 

5 17 44.43 

163 01 29.88 

7 

C-5 

-8 05 40.58 

156 49 24.82 

49 1 

Provisional 

DOS 

-8 05 40.58 

156 49 24.82 

49 

C-5 

1 21 42.13 

172 55 47.26 

7* 

1966 SECOR 
ASTRO 

1 21 42.13 

172 55 47.26 

7 

C-5 

-17 45 31.01 

177 27 02.83 

17 f 

Viti 

Levu 1916 

-17 45 31.01 

177 27 02.83 

17 

C-5 

-2 46 28.99 

188 16 43.47 

6t 

1966 Canton 
ASTRO 

-2 46 28.99 

188 16 43.47 

6 

C-5 

16 43 51.68 

190 28 41.55 

6t 

Johnston 
Island 1961 

16 43 51.68 

190 28 41.55 

6 

C-5 

28 12 32.06 

182 37 49.53 

6 

Midway 
ASTRO 1961 

28 12 32.06 

182 37 49.53 

6 

C-5 

20 49 37.00 

203 31 52.77 

32 

Old Hawaiian 

20 49 25.14 

203 32 01.88 

31 

C-5 

31 55 18.41 

278 26 00.26 

29 

N.A. 

31 55 18.86 

278 26 00.18 

-27 

C-5 

32 00 04.04 

278 50 43.17 

27 

NA. 

32 00 04.49 

278 50 43.13 

-30 

C-5 

25 29 21.18 

279 37 39.35 

18 

NA. 

25 29 22.51 

279 37 39.37 

-44 

. 

C-5 


51 



Table D8 



BELTVL 6002 39°01' 39703 283‘T0’26'. , 94 45 NA. 

39 01 38.80 283 10 27.40 -10 C-5 

ASTRMD 6100 39 01 39.72 283 10 27.83 45 N A. 

39 01 39.49 283 10 28.29 -10 C-5 

TIM INS 6113 48 33 56.17 278 37 44.54 290 NJV. 

48 33 55.70 278 37 44.49 232 C-5 


Table D9 

SPEOPT— Optical— Source B. 


Name 

Station No. 

Latitude 

Longitude 

Geodetic 

Height 

(meters) 

Datum 

1UNDAK 

7034 

48°01'21 T . r 40 

262°59 T 2rj56 

255 

N.A. 



48 01 20.81 

262 59 19.55 

201 

C-5 

1EDINB 

7036 

26 22 45.44 

261 40 09.03 

67 

N.A. 



26 22 46.35 

261 40 07.34 

15 

C-5 

1COLBA 

7037 

38 53 36.07 

267 47 42.12 

271 

N.A. 



38 53 35.81 

267 47 40.85 

218 

C-5 

1BERMD 

7039 

32 21 48.83 

295 20 32.56 

21 

NA. 



32 21 48.94 

295 20 34.18 

-28 

C-5 

1PURIO 

7040 

18 15 26.22 

294 00 22.17 

58 

N.A. 



18 15 28.30 

294 00 23.63 

5 

C-5 

1GSFCP 

7043 

39 01 15.01 

283 10 19.93 

54 

NA. 



39 01 14.78 

283 10 20.39 

-1 

C-5 

1CKVLE 

7044 

38 22 12.50 

274 21 16.81 

187 

N.A. 



38 22 12.33 

274 21 16.28 

131 

C-5 

1DENVR 

7045 

39 28 48.03 

255 23 41.19 

1796 

N .A. 



39 38 47.54 

255 23 38.52 

1751 

C-5 

1JUM24 

7071 

27 01 12.77 

279 53 12.31 

25 

N.A. 



27 01 14.00 

279 53 12.30 

-38 

C-5 

1JUM40 

7072 

27 01 13.17 

279 53 12.49 

25 

N.A. 



27 01 14.39 

279 53 12.49 

-38 

C-5 

1JUPC1 

7073 

27 01 13.11 

279 53 12.72 

22 

N.A. 



27 01 14.33 

279 53 12.72 

-41 

C-5 

1JUBC4 

7074 

27 01 13.33 

279 53 12.76 

25 

N.A. 



27 01 14.55 

279 53 12.76 

-38 

C-5 

1SUDBR 

7075 

46 27 20.99 

279 03 10.35 

281 

NA. 



46 27 20.52 

279 03 10.35 

224 

C-5 

1JAMAC 

7076 

18 04 31.98 

283 11 26.52 

485 

N.A. 



18 04 34.20 

283 11 27.03 

423 

C-5 



Table DIO 

SPEOPT— Laser— Source B 


Name 

Station No. 

Latitude 

GODLAS 

7050 

39°01 , 13'. l 68 



39 01 13.45 

ROSLAS 

7051 

35 11 46.60 



35 11 46.70 


Longitude 

Geodetic 

Height 

(meters) 

Datum 

assno'isvos 

55 

N A. 

283 10 18.51 

0 

C-5 

277 07 26.23 

879 

NA. 

277 07 26.02 

822 

C-5 


Table Dll 

International— Optical— Source I.* 







Geodetic 


Source 

Name 

Station No. 

Latitude 

Longitude 

Height 

Datum 






(meters) 



DELFTH 

8009 

52°00 , 09’. , 24 

4°22 , 21"23 

23 

European 




52 00 06.12 

4 22 15.30 

28 

C-5 

D 

ZIMWLD 

8010 

46 52 41.77 

7 27 57.56 

903 

BERNE 




46 52 36.73 

7 27 52.54 

907 

C-5 


MALVRN 

8011 

52 08 39.12 

358 01 59.49 

111 

European 




52 08 35.68 

358 01 53.03 

125 

C-5 


*UnIess "Source” is specified otherwise. 


Table D12 
SAO— Optical. 


Name 

Station No. 

Location 

lORGAN 

9001 

Organ Pass, New Mexico 

lOLFAN 

9002 

Olifantsfontein, South Africa 

WOOMER 

9003 

Woomera, Australia 

1 SPAIN 

9004 

San Fernando, Spain 

1TOKYO 

9005 

Tokyo, Japan 

1NATOL 

9006 

Naini Tal, India 

1QUIPA 

9007 

Arequipa, Peru 

1SHRAZ 

9008 

Shiraz, Iran 

1CURAC 

9009 

Curacao, Lesser Antilles 

1JUPTR 

9010 

Jupiter, Florida 

1VILDO 

9011 

Villa Dolores, Argentina 

1MAUIO 

9012 

Maui, Hawaii 

OSLONR 

9426 

Oslo, Norway 

AUSBAK 

9023 

Woomera, Australia 

NATALB 

9029 

Natal, Brazil 

AGASSI 

9050 

Cambridge, Massachusetts 

COLDLK 

9424 

Cold Lake, Alberta 

EDWAFB 

9425 

Edwards AFB, California 

RIGLAT 

9428 

Riga, Latvia 

POTDAM 

9429 

Potsdam, Germany 

ZVENIG 

9430 

Zvenigorod, Russia 
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Table D13 
STADAN— Optical 


Name 

Station No. 

Location 

1BPOIN 

1021 

Blossom Point, Maryland 

1FTMYR 

1022 

Fort Myers, Florida 

lOOMER 

1024 

Woomera, Australia 

lQurro 

1025 

Quito, Ecuador 

1LIMAP 

1026 

Lima, Peru 

ISA TAG 

1028 

Santiago, Chile 

1MOJAV 

1030 

Mojave, California 

1JOBUR 

1031 

Johannesburg, Union of South Africa 

1NEWFL 

1032 

St. John’s, Newfoundland 

1COLEG 

1033 

College, Alaska 

1GFORK 

1034 

East Grand Forks, Minnesota 

1WNKFL 

1035 

Winkfield, England 

1ROSMA 

1042 

Rosman, North Carolina 

1TANAN 

1043 

Tananarive, Madagascar 


_ 
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Table D14 
STADAN— R/R. 


Name 

Station No. 

Location 

MADGAR 

ROSRAN 

CARVON 

1122 

1126 

1152 

Tananarive, Madagascar 
Rosman, North Carolina 
Carnarvon, Australia 

Table D15 

Navy T R A NE T — -Dopple r . 

Name 

Station No. 

Location 

LASHAM 

2006 

Lasham, England 

SANHES 

2008 

Sao Jose dos Campos, Brazil 

PHILIP 

2011 

San Miquel, Philippines 

SMTHFD 

2012 

Smithfield, Australia 

MISAWA 

2013 

Misawa, Japan 

ANCHOR 

2014 

Anchorage, Alaska 

TAFUNA 

2017 

Tafuna, American Samoa 

THULEG 

2018 

Thule, Greenland 

MCMRDO 

2019 

McMurdo Sound, Antarctica 

WAHIWA 

2100 

South Point, Hawaii 

LACRES 

2103 

Las Cruces, New Mexico 

LASHM2 

2106 

Lasham, England 

APLMND 

2111 

APL Howard County, Maryland 

PRETOR 

2115 

Pretoria, Union of South Africa 

SHEMYA 

2739 

Shemya Island, Alaska 

BELTSV 

2742 

Beltsville , Maryland 

STNVIL 

2745 

j Stoneville, Mississippi 
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Table D16 


Air Force— Optical. 


Name 

Station No. 

Location 

ANTIGA 

3106 

Antigua Island, Lesser Antilles 

GRNVLE 

3333 

Stoneville, Mississippi 

GRVILL 

3334 

Stoneville, Mississippi 

USAFAC 

3400 

Colorado Springs, Colorado 

BEDFRD 

3401 

L. G. Hanscom Field, Massachusetts 

SEMMES 

3402 

Semmes Island, Georgia 

SWANIS 

3404 

Swan Island, Caribbean Sea 

GRDTRK 

3405 

Grand Turk, Caicos Islands 

CURACO 

3406 

Curacao, Lesser Antilles 

TRNDAD 

3407 

Trinidad Island 

GRANFK 

3451 

Grand Forks, North Dakota 

TWINOK 

3452 

Twin Oaks, Oklahoma 

ROTHGR 

3453 

Rothwesten, West Germany 

ATHNGR 

3463 

Athens, Greece 

TORRSP 

3464 

Torre j on de Ardoz, Spain 

CHOFUJ 

3465 

Chofu, Japan 

KINDLY 

3471 

Kindley AFB, Bermuda 

HUNTER 

3648 

Hunter AFB, Georgia 

JUPRAF 

3649 

Jupiter, Florida 

ABERDN 

3657 

Aberdeen, Maryland 

HOMEST 

3861 

Homestead AFB, Florida 

CHYWYN 

3902 

Cheyenne, Wyoming 


Table D17 

Army Map Service — SECOR. 


Name 

Station No. 

Location 

HERNDN 

5001 

Herndon, Virginia 

CUBCAL 

5200 

San Diego, California 

LARSON 

5201 

Moses Lake, Washington 

WRGTON 

5202 

Worthington, Minnesota 

GREENV 

5333 

Greenville, Mississippi 

TRUKIS 

5401 

Truk Island, Caroline Islands 

SWALLO 

5402 

Swallow Island, Santa Cruz Islands 

KUSAIE 

5403 

Kusaie Island, Caroline Islands 

GIZZOO 

5404 

Gizzoo, Gonzongo, Solomon Islands 

TARAWA 

5405 

Tarawa, Gilbert Islands 

NANDIS 

5406 

Nandi, Viti Levu, Fiji Islands 

CANTON 

5407 

Canton Island, Phoenix Islands 

JONSTN 

5408 

Johnston Island, Pacific Ocean 

MIDWAY 

5410 

Eastern Island, Midway Islands 

MAUIHI 

5411 

Maui, Hawaii 

FTWART 

5648 

Fort Stewart, Georgia 

HNTAFB 

5649 

Hunter AFB, Georgia 

HOMEFL 

5861 

Homestead AFB, Florida 
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Table D18 


USC &GS — Optical . 


Name 

Station No. 

Location 

BELTVL 

6002 

Beltsville, Maryland 

ASTRMD 

6100 

Beltsville, Maryland 

TIM INS 

6113 

Timmins, Ontario 


Table D19 
S PE OPT— Optical. 


Name 

Station No. 

Location 

1UNDAK 

7034 

Univ. North Dakota, Grand Forks, 
North Dakota 

1EDINB 

7036 

Edinburg, Texas 

1COLBA 

7037 

Columbia, Missouri 

1BERMD 

7039 

Bermuda Island 

1PURIO 

7040 

San Juan, Puerto Rico 

1GSFCP 

7043 

GSFC, Greenbelt, Maryland 

1CKVLE 

7044 

Clarksville, Indiana 

1DENVR 

7045 

Denver, Colorado 

1JUM24 

7071 

Jupiter, Florida 

1JUM40 

7072 

Jupiter, Florida 

1JUPC1 

7073 

Jupiter, Florida 

1JUBC4 

7074 

Jupiter, Florida 

1SUDBR 

7075 

Sudbury, Ontario 

1 JAM AC 

7076 

Jamaica, B. W. I. 


Table D20 
SPEOPT— Laser. 


Name 


GODLAS 

ROSLAS 


Station No. 

7050 

7051 


Location 

GSFC, Greenbelt, Maryland 
Rosman, North Carolina 


Table D21 


International —Optical . 


Name 


DELFTH 

ZIMWLD 

MALVRN 


Station No. 

8009 

8010 
8011 


Location 
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Delft, Holland 
Berne, Switzerland 
Malvern, England 
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Appendix E 

Table of GEOS I Orbit Numbers from 
12/31/65 to 1/5/66 
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Table El 


Pass Numbers for GEOS I. 

(Mean Period of GEOS I at epoch 12/31/65 01 H 38 M 22 s = 29005393 = 12093236) 


Pass No. 

From 

To 

650 

12/31/65 

01 H 

18 m 

12/31/65 

03 H 

17 M 

651 


03 

18 


05 

17 

652 


05 

18 


07 

18 

653 


07 

19 


09 

18 

654 


09 

19 


11 

18 

655 


11 

19 


13 

18 

656 


13 

19 


15 

19 

657 


15 

20 


17 

19 

658 


17 

20 


19 

19 

659 


19 

20 


21 

20 

660 


21 

21 


23 

20 

661 


23 

21 

01/01/66 

01 

20 

662 

01/01/66 

01 

21 


03 

21 

663 


03 

22 


05 

21 

664 


05 

22 


07 

21 

665 


07 

22 


09 

21 

666 


09 

22 


11 

22 

667 


11 

23 


13 

22 

668 


13 

23 


15 

22 

669 


15 

23 


17 

23 

670 


17 

24 


19 

23 

671 


19 

24 


21 

23 

672 


21 

24 


23 

24 

673 


23 

25 

01/02/66 

01 

24 

674 

01/02/66 

01 

25 


03 

24 

675 


03 

25 


05 

24 

676 


05 

25 


07 

25 

677 


07 

26 


09 

25 

678 


09 

26 


11 

26 

679 


11 

26 


13 

26 

680 


13 

27 


15 

26 

681 


i ~ 

27 


17 

26 

682 


17 

27 


19 

27 

683 


19 

28 


21 

27 

684 


21 

28 


23 

27 

685 


23 

28 

01/03/66 

01 

28 

686 

01/03/66 

01 

29 


03 

28 

687 


03 

29 


05 

28 

688 


05 

29 


07 

28 

689 


07 

29 


09 

29 

690 


09 

30 


11 

29 

691 


11 

30 


13 

29 

692 


13 

30 


15 

30 

693 


15 

31 


17 

30 

694 


17 

31 


19 

30 

695 


19 

31 


21 

31 

696 


21 

32 


23 

31 

697 


23 

32 

01/04/66 

01 

31 


61 



\ 


Table El— Continued 


Pass No. 


From 



To 


698 

0 1/04/66 

01 H 

32 M 

01/04/66 

03 H 

31 M 

699 


03 

32 


05 

32 

700 


05 

33 


07 

32 

701 


07 

33 


09 

32 

702 


09 

33 


11 

33 

703 


11 

34 


13 

33 

704 


13 

34 


15 

33 

705 


15 

34 


17 

34 

706 


17 

35 


19 

34 

707 


19 

35 


21 

34 

708 


21 

35 


23 

34 

709 


23 

35 

01/05/66 

01 

35 

710 

01/05/66 

01 

36 


03 

35 

711 


03 

36 


05 

35 

712 


05 

36 


07 

36 
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